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Abstract 


We perform a series of stationary wave model (SWM) experiments in which the boreal summer 
atmosphere is forced, over a number of locations in the continental U.S., with an idealized 
diabatic heating anomaly that mimics the atmospheric heating associated with a dry land surface. 
For localized heating within a large portion of the continental interior, regardless of the specific 
location of this heating, the spatial pattern of the forced atmospheric circulation anomaly (in 
terms of 250-mb eddy streamfunction) is largely the same: a high anomaly forms over west- 
central North America and a low anomaly forms to the east. In supplemental atmospheric 
general circulation model (AGCM) experiments, we find similar results; imposing soil moisture 
dryness in the AGCM in different locations within the US interior tends to produce the 
aforementioned pattern, along with an associated near-surface warming and precipitation deficit 
in the center of the continent. The SWM-based and AGCM-based patterns generally agree with 
composites generated using reanalysis and precipitation gauge data. The AGCM experiments 
also suggest that dry anomalies imposed in the lower Mississippi Valley have remote surface 
impacts of particularly large spatial extent, and a region along the eastern half of the US-Canada 
border is particularly sensitive to dry anomalies in a number of remote areas. Overall, the SWM 
and AGCM experiments support the idea of a positive feedback loop operating over the 
continent: dry surface conditions in many interior locations lead to changes in atmospheric 


circulation that act to enhance further the overall dryness of the continental interior. 
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1. Introduction 


From a societal perspective, much of what characterizes the Earth’s climate can be 
related to atmospheric processes — the mean and variability of rainfall in a region, for example, 
or the statistics of near-surface air temperature. The atmospheric component of climate, with its 
winds, storms, clouds, and aerosols, along with myriads of other processes, is overwhelmingly 
complex in and of itself. The full climate system, however, is made even more complex by the 
interactions of the atmosphere with the system’s other components, notably the ocean, the land, 
and the cryosphere. The coupling of these different components allows, for example, the 
development and maintenance of modes of behavior, such as the El Nino / Southern Oscillation 


(ENSO) phenomenon, that otherwise would not exist. 


While the coupling of the different climate components increases the overall complexity 
of the system, it also has the beneficial impact of sometimes increasing climate predictability — 
the ability to predict ahead of time, perhaps weeks to seasons in advance, a climatic variation 
such as a regional drought. Atmospheric processes are strongly influenced by chaotic 
atmospheric dynamics. Because of chaos in the atmosphere, forecasts relying solely on 
atmospheric initialization cannot expect to have skill exceeding about two weeks. Time scales of 
variation in the ocean, land, and sea ice, however, and of various coupled modes such as ENSO, 
are much longer, and these longer timescales can translate, through coupling, to predictive skill 
for various atmospheric quantities. Operational seasonal forecasting systems indeed rely on the 
added predictability associated with the coupling of the atmosphere to slower components of the 


system (NRC, 2010). 


The coupling of the land to the atmosphere potentially allows the long time scales of soil 


moisture anomalies (weeks to months; Entin et al. 2000; Vinnikov and Yeserkepova 1991) and 
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snow anomalies (winter through the spring melt season) to contribute to atmospheric 
predictability at these time scales. The idea is simple — if a soil moisture anomaly, for example, 
is known at the start of a forecast, it can be predicted with some skill weeks to months into the 
forecast due to the slow time scale over which it varies. If, weeks to months into the forecast, the 
atmosphere in turn responds in a predictable way to the predicted soil moisture anomaly, then 
some skill is imparted to atmospheric prediction at these leads. The impacts of soil moisture 
variations on atmospheric variability have been studied extensively using both climate model 
analyses (e.g., Delworth and Manabe 1989, Douville and Chauvin 2000, Koster et al. 2006, Guo 
et al. 2006, Dirmeyer et al. 2013) and observational analyses (e.g., Betts and Ball 1995, Findell 
and Eltahir 1997, Koster et al. 2003, Taylor et al. 2011). Positive impacts of realistic soil 
moisture initialization on precipitation and air temperature forecast skill at monthly leads were 
quantified across a broad range of forecasting systems in the Global Land-Atmosphere Coupling 
Experiment, Phase 2 (Koster et al. 2011). An impact of snow initialization on forecast skill was 


demonstrated by Peings et al. (2011) and Ambadan et al. (2015). 


The present paper focuses specifically on soil moisture anomalies, their impacts on the 
overlying atmosphere, and their associated impacts on remote precipitation and air temperature 
fields. The physical mechanisms underlying these impacts are not fully known. A local effect 
on the air temperature is straightforward; higher soil moisture levels can lead to increased 
evapotranspiration, which in turn leads to increased evaporative cooling of the land surface and, 
consequently, of the overlying air (Seneviratne et al. 2010). A local impact on precipitation can 
also be envisioned: higher soil moistures can increase evapotranspiration and accordingly reduce 
surface sensible heat flux, leading to a modification of the planetary boundary layer and thus to a 


modification of conditions determining the onset of moist convection (Betts et al. 1994). To 
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what extent, however, can a soil moisture anomaly affect meteorological conditions in remote 
locations? On this subject the literature is more limited. For such a remote impact to be 
predictable — for information at the site of the soil moisture anomaly to be translated to the 
remote region — atmospheric transports or the atmospheric circulation itself would need to be 
modified in a predictable way. Douville et al. (2002) suggested that certain soil moisture 
anomalies may help trigger stationary planetary waves over Europe that in turn may affect the 
transport of dry air from midlatitudes to the tropics. Lau and Kim (2012) connected extreme 
weather events in Russia and Pakistan through atmospheric patterns that were potentially 
amplified by dry land conditions. Koster et al. (2014) showed that forcing an atmospheric 
general circulation model (AGCM) with a soil moisture dipole over the continental US — wet 
conditions in the Pacific Northwest and dry conditions in the southern Great Plains — leads to a 
distinct planetary wave pattern that induces further drying and warming in the latter region. 
Taylor et al. (2011) demonstrate the impact of surface moisture heterogeneity on the formation 
of mesoscale circulations that promote storm formation over dry regions. It is worth noting that 
the impact of other land surface anomalies on the atmospheric circulation and associated remote 
impacts have also been examined; both Robock et al. (2003) and Cohen et al. (2014), for 
example, examined the role of Siberian snow cover on the atmospheric circulation and the 
associated potential for prediction, and Xue et al. (2012) examined the ability of subsurface soil 
temperature anomalies in the western US to affect conditions in the eastern US via Rossby 


waves. 


In the present paper we expand on the findings of Koster et al. (2014) through a more 
comprehensive analysis of the impacts of North American soil moisture anomalies on the 


overlying atmospheric circulation and remote meteorological variables, with a focus on boreal 
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summer. Results are shown in the order of the complexity of the models used to produce them. 
That is, after describing the models and the datasets examined (section 2), we first show a set of 
results (section 3a) obtained by forcing a relatively simple stationary wave model (SWM) with 
an idealized heating anomaly that mimics the direct effect of surface drying. The SWM results 
have the advantage of demonstrating clearly one of our main findings, namely, that the 
atmospheric circulation responds in a similar way to heating anomalies imposed in very different 
locations. The SWM results, of course, also have the disadvantage of lacking physical processes 
and hence the physical complexities of a full AGCM; our second set of results (section 3b) are 
therefore obtained with a full suite of corresponding AGCM experiments. The AGCM results 
support the aforementioned main finding, though with the expected increase in noise and with 
added nuances in the patterns produced. The AGCM experiments also provide information on 
remote precipitation and air temperature impacts. We then examine available observations-based 
data (section 3c) and show that they in turn provide support for the SWM and AGCM findings. 
We conclude (section 3d) with supplemental analyses of the AGCM data that illustrate how the 


strength of soil moisture-meteorology teleconnections may vary in space. 


Our AGCM analysis, by the way, utilizes large ensembles (~100-200 members), a 
reflection of the subtlety of the signals we seek here. We are, in effect, seeking the net impact 
(e.g., on mean states) of shifts in the probability density functions of various quantities rather 
than more deterministic relationships that would necessarily lead to a first-order improvement in 


predictions. This caveat should be kept in mind when evaluating the analyses that follow. 


2. Models and Datasets Used 
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a. Stationary Wave Model (SWM) 


The nonlinear and time-dependent stationary wave model (SWM) used in this study is 
based on the three-dimensional primitive equations in o coordinates. All of the basic variables in 
the model are deviations from a prescribed zonal mean or 3-D climatological flow. The basic 
prognostic equations are those for perturbation vorticity, divergence, temperature, and the 
logarithm of surface pressure. The model-generated transient disturbances are suppressed by 
applying strong damping. The stationary wave solution is obtained by integrating the model to a 


quasi-steady state after a short period of time. 


The model has rhomboidal wavenumber-30 (R30) truncation in the horizontal and 14 
unevenly spaced o levels in the vertical. While the R30 resolution (roughly 2.25° latitude x 
3.75° longitude) is coarser than that of the AGCM described below (roughly 1° resolution), we 
nevertheless expect the AGCM and SWM results to be directly comparable. This is because, 
relative to AGCM results, SWM results are not as affected by resolution. Much of the 
resolution-dependent behavior in an AGCM stems from its use of model physical 
parameterizations to produce stationary wave forcings such as diabatic heating. In contrast, 
diabatic heating in the SWM is prescribed and is thus not resolution-dependent. Output from our 
R30 SWM experiments and that from SWM experiments at roughly 1° resolution should be 


similar. 


The SWM has been shown to be a valuable tool for diagnosing the maintenance of the 


climatological atmospheric circulation and its anomalies (in terms of magnitude and pattern) as 
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well as for investigating dynamics of stationary waves (e.g., Ting et al. 2001; Held et al. 2002; 


Schubert et al. 2011). See Ting and Yu (1998) for more details on the SWM. 


b. Atmospheric General Circulation Model (AGCM) 


The simulations examined here were performed with the GEOS-5 (Global Earth 
Observing System, Version 5) atmospheric general circulation model (AGCM), a component of 
the state-of-the-art modeling system maintained by the Global Modeling and Assimilation Office 
(GMAO) of the National Aeronautics and Space Administration’s Goddard Space Flight Center 
(NASA/GSFC). The version of the GEOS-5 AGCM used here is run on a cubed-sphere grid 
(Putman and Lin 2007). As is standard for this kind of system, unresolvable physical processes 
are parameterized; GEOS-5 uses the approaches of Moorthi and Suarez (1992) for convection, 
Bacmeister et al. (2006) for prognostic cloud cover, Chou and Suarez (1994) and Chou 
(1990,1992) for longwave and shortwave radiation processes, Lock et al. (2000) for turbulence 
near the surface, and Helfand and Schubert (1995) for the surface layer (though over land the 
surface layer includes a viscous sublayer). Land surface processes are simulated with the 
Catchment land surface model (Koster et al. 2000), a model that simulates explicitly the impact 


of topography-driven subgrid soil moisture variability on the surface energy and water balances. 


As with all AGCMs, the GEOS-5 AGCM has biases in its simulated climate. 
Nevertheless, the AGCM is well-vetted. For the past several years, short-term weather forecasts 
produced by the model have been analyzed extensively by GMAO personnel every week to 
maintain confidence in model performance and to determine and correct deficiencies; this 


process ensures the continued maintenance of a state-of-the-art system. The AGCM has been 
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deemed suitable for use in GMAO operations, serving, for example, as the model component of 
the MERRA and MERRA-2? reanalyses (Modern-Era Retrospective Analysis for Research and 
Applications; Rienecker et al. 2011; Bosilovich et al. 2016) and, when coupled to an ocean 
model, as the basis for GMAO seasonal forecasts (Ham et al. 2014). The GEOS-5 AGCM has 
been used extensively to study the mechanisms underlying climate variability (Schubert et al. 
2014, Koster et al. 2015, Wang et al. 2016). Molod et al. (2012) document the performance 
(including biases) of the free-running model, and Gelaro et al. (2014) and Bosilovich et al. 
(2016) show, respectively, how the model behaves at very high resolution and within an analysis 


system. 


c. Observations-Based Data 


The SWM and AGCM simulations produce streamfunction fields, and the AGCM also 
produces near-surface air temperatures and precipitation rates. To evaluate the realism of these 


model results, we interpret them in the context of available observations-based data. 


We use reanalysis data to provide observations-based information on soil moisture, 
streamfunction, and near-surface air temperature. In the MERRA-2 reanalysis (Bosilovich et al. 
2016), the land surface is forced with precipitation rates that are corrected with gauge-based 
precipitation observations (Reichle and Liu 2014). As a result, soil moisture contents in the 
MERRA-? reanalysis evolve realistically, attaining the level of accuracy seen in traditional Land 
Data Assimilation System, or LDAS, operational products (e.g., Xia et al 2014). Note that 
because we utilize MERRA-2 root zone soil moistures only to determine periods for which the 


real world experienced dry land surface conditions, their level of accuracy should be acceptable; 
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while the products of different LDAS systems do show differences, they nevertheless tend to be 
highly consistent in their identification of extreme dry periods (e.g., Koster et al. 2009, Wang et 


al. 2009). 


We could also use MERRA-? to provide the observations-based streamfunction and air 
temperature data we need, but we choose instead to extract these data from the European Centre 
for Medium-Range Weather Forecasts (ECMWF) ERA-Interim reanalysis (Dee et al. 2011). 
This choice is driven by the fact that the AGCM experiments described below and the MERRA- 
2 reanalysis utilize the same atmospheric model, which could potentially promote similar biases 
in both sets of data. ERA-Interim provides observations-based data that are inherently more 
independent. (As it turns out, the use of MERRA-? produces essentially the same results [not 
shown].) In addition, the near-surface air temperatures produced by ERA-Interim are known to 
be of high quality (Simmons et al. 2010), probably in part because ERA-Interim uses station 


measurements of air temperature to update soil temperatures. 


For precipitation, we use the observations-corrected precipitation output from MERRA-?2. 
Because we only look at precipitation totals for 10-day periods in June and July, the use of the 
observations-corrected MERRA-?2 data is equivalent to using, over land, precipitation data from 
the Climate Prediction Center (CPC) Unified Gauge-Based Analysis of Global Daily 
Precipitation (CPCU; ftp://ftp.cpc.ncep.noaa.gov/precip/CPC_UNI_ PRCP/GAUGE _GLB/) 
database and, over the ocean, from a mix of data from the Global Precipitation Climatology 
Project (GPCP, Adler et al. 2003, Huffman et al. 2009) and the CPC Merged Analysis of 


Precipitation (CMAP; ftp://ftp.cpc.ncep.noaa.gov/precip/cmap/). In essence, then, we compare 


our AGCM precipitation results over land directly to rain gauge measurements. See Reichle and 
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Liu (2014) for a description of how the different observational precipitation datasets were 


incorporated into the MERRA-2 data. 


3. Experiments and Results 


a. SWM Responses to Idealized Diabatic Heating Anomalies 


In nature, a drier-than-average surface state (i.e., low soil moisture) can affect diabatic 
heating in the overlying atmosphere through its impact on the surface energy balance — the dry 
conditions lead to an anomalously high sensible heat flux from the surface and thus to increased 
near-surface heating in the immediate area of the anomaly. When, in contrast, the surface is 
anomalously wet and latent heat flux dominates, the associated anomalous diabatic heating may 
appear geographically far away, depending on where that excess evaporated moisture eventually 


condenses. 


To represent the immediate atmospheric impact of a dry surface anomaly and the 
associated local increase in sensible heat flux, we impose in the SWM, over a selected area of 
dryness, an idealized diabatic heating anomaly with the vertical profile shown in Figure la. The 
areal extent of the imposed diabatic heating anomaly is given horizontal half-widths of 5° 
longitude and 5° latitude, as indicated in Figure 1b for one of the experiments. This vertical 
distribution and spatial extent is consistent with the structures of the diabatic heating anomalies 
produced in the AGCM experiments described in section 3b below. To illustrate this, the 
vertical and spatial (zonal) structure of diabatic heating produced in a representative AGCM 
experiment (specifically, AGCM-L, to be defined below) is illustrated in Figure Ic. 
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In a standard experiment, the anomalous heating is imposed throughout the model 
integration. The mean basic state used for the SWM is the 3-D June-July averaged climatological 
(1980-2010) basic flow taken from the GEOS-5 AGCM. The SWM reaches a steady state at 
about day 25, meaning that results from any day thereafter would represent the steady-state 


model response. Here we present the model eddy streamfunction solution at day 50. 


We use this approach in a series of 21 experiments with the SWM, with each experiment 
assuming a different assumed area of dryness and thus a different geographical placement of the 
idealized anomaly shown in Figure 1. Figure 2 shows the 21 geographical areas considered. Ina 
given experiment, the imposed idealized anomalous diabatic heating was centered as close as 
possible (given the grid resolution of the SWM) to the central longitude and latitude of that 
experiment’s indicated area. The 21 SWM experiments are hereafter labeled SWM-A through 


SWM-U. 


Figure 3 shows the response of the 250-mb eddy streamfunction field to the imposed 
heating anomaly in each experiment. Notice the similarity in the responses. In almost all of the 
experiments, a high anomaly is seen in the west-central part of the northern half of the continent, 
and in most of the experiments, a low anomaly is seen further to the east. The eddy 
streamfunction patterns generally look very similar in their placement, though the magnitudes of 
the responses are substantially reduced when the diabatic heat source is placed on the western 
coast, and the patterns do change when it is placed toward the eastern coast. Certainly within the 
continental interior, different placements of the heat source promote a common atmospheric 
circulation response, with the maximum and minimum of the streamfunction field not simply 


shifting laterally with the heat source. 
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We choose the following approach to characterize the degree of similarity amongst the 
responses. We first average the 21 streamfunction anomaly fields in Figure 3 into a single field 
and then compute the square of the spatial correlation, r’, between each experiment’s 250-mb 
streamfunction anomaly field and this average field over the North American area shown. The 
heavy dashed lines in Figure 4a enclose the nine experiments for which this r’ is at least 0.65, 
i.e., for which the average field over the 21 experiments explains about 2/3 or more of the spatial 
variance of a given experiment’s field. Figure 4b then presents the results of a second averaging 
exercise; it shows the results of averaging the 250-mb streamfunction anomaly fields over the 


nine experiments indicated in Figure 4a. 


A comparison of the individual streamfunction anomaly fields in Figure 3 with the 
pattern in Figure 4b shows that this pattern is indeed a first-order representation of the SWM 
results for these nine experiments. This point is underlined by Figure 4c, which shows the 
square of the correlation coefficient between each experiment’s streamfunction anomaly field 
and the 9-experiment average in Figure 4b. Each of the 9 experiments focusing on heating 
within the delimited area produces a streamfunction field in strong agreement with the 9- 


experiment average, with r* values upward of 0.7. 


There are presumably many ways to characterize similarity in experiment results. Based 
on the approach used here, we make the claim that the response of the 250-mb streamfunction 
field to diabatic heating anywhere within the region delimited in Figure 4a is similar in pattern — 
regardless of where surface heating is located in this region, the atmospheric circulation responds 


in much the same way. Much of our later analyses will focus on this interior region. 
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b. AGCM Responses to Imposed Land Surface Dryness 


The next step in our analysis is to determine if the basic pattern (Figure 4b) found by the 
SWM for the interior region in Figure 4a is also found within a full suite of AGCM experiments. 


We describe the AGCM experiments below. 


(i) Control Ensemble. The control for our experiments consists of an ensemble of 768 4-month 
simulations with the atmosphere and land components of GEOS-5. This decidedly large number 
of ensemble members helps ensure accurate statistics. Initialization followed the strategy 
employed by Koster et al. (2014): each simulation was initialized on 1 April 2012, with initial 
atmosphere and land conditions taken from different years of the MERRA reanalysis, and with 
the initial atmospheric conditions perturbed slightly to produce multiple ensemble members for 
each MERRA year. Sea surface temperatures for 2012 (a year for which the real world 
experienced warming in the Great Plains of the US) were prescribed using Atmospheric Model 
Intercomparison Study (AMIP)-style protocols (Gates 1992). Each simulation was performed on 
a C90 cubed-sphere grid (equivalent to 1°x1° resolution), and output from each simulation was 


subsequently written out onto a 1°x1° latitude-longitude grid for analysis. 


(ii) Experiment Ensembles. A total of 21 experiments were performed, each experiment 
consisting of an ensemble of either 96 or 192 simulations that differed from the control only in 
the prescribed drying of a selected land area. The drying of the selected area proceeded as 
follows. During every time step from | April to 30 June, any precipitation produced over the 


selected area by the AGCM was artificially zeroed before it reached the land surface, thereby 
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allowing soil moistures in the land model there to evolve to levels corresponding to a 3-month 
dry period. The artificial zeroing of the precipitation ceased on 1 July, so that during the month 
of July, soil moistures were allowed to recover. The zeroing of the precipitation during 1 April - 
30 June was only at the surface; the AGCM’s atmospheric vapor and temperature fields, having 
been modified by condensation processes in the production of precipitation, were not artificially 


reset to their pre-rainfall values. 


The 21 areas considered in the AGCM experiments mimic those used for the SWM 
experiments (see Figure 2), and we will thus refer to the AGCM experiments as AGCM-A 
through AGCM-U. The regions examined with the AGCM are, in fact, slightly smaller, an 
artifact of mapping issues associated with our use of a “cube grid” for the AGCM, a grid that 
does not translate exactly to a latitude-longitude grid. The areas, each roughly of size 7°x7°, 
essentially span the continental US. While the individual areas are somewhat smaller than the 
two examined by Koster et al. (2014), dryness anomalies within them are still able to affect the 
AGCM’s overlying circulation, as will be shown later. The number of ensemble members 


comprising each experiment is listed in parentheses below the experiment identifier in Figure 2. 


As in the SWM analysis, experiment results relative to the control are analyzed for the 
period 1 June - July 31, a time when land surface impacts on the atmosphere are expected to be 
large due to high surface turbulent fluxes. (We indeed expect land impacts to be reduced during 
fall through spring, when the turbulent fluxes are reduced due to reduced solar forcing.) In the 
context of the experimental design, the period examined corresponds to the final month of an 
extreme 3-month meteorological drought (June) and the first recovery month thereafter (July). 
Note that while a precipitation deficit was not imposed in the specified area during July, 


persistence ensured that soil moisture levels there during July remained significantly low. This 
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is illustrated in Figure 5a, which shows the percentile, based on the statistics of the control 
ensemble, of average July root zone soil moisture content produced by the experiments at each 
grid cell. The map is in fact a composite of results from the different experiments; percentiles 
shown within a given area are from the experiment that utilized that area. (High values along 
some edges, by the way, reflect the aforementioned slight inconsistency between the cube grid 
and a regular latitude-longitude grid.) July soil moistures tend to lie below the 5" percentile, 
indicating that soil moistures in the prescribed dryness regions, forced to be dry in June, do 
indeed remain dry through July. Deserts are an obvious exception, with the percentiles within 
Region P in particular sometimes strongly (and surprisingly) exceeding 50%; note, however, that 
these regions feature a very small range of soil moisture in the control, with all control values at 
the very dry end. The high percentiles in Region P are presumably artifacts, and the high values 


notwithstanding, July soil moistures in Region P for Experiment P are themselves very low. 


This last point is emphasized in Figure 5b, which shows the difference in the total 
evaporation for June and July between each experiment and the control. As in Figure 5a, the 
map in Figure 5b is a composite; results plotted within a given region are for the experiment 
focusing on that region (and indeed are derived, for this plot, from the first 96 ensemble 
members of the experiment). July evaporation rates from Region P in experiment AGCM-P are 
not significantly increased, despite the percentiles shown for the region in Figure 5a. More 
broadly, Figure 5b is a useful reference because it indicates the relative strengths of the local 
diabatic heating anomaly effectively imposed in the different AGCM experiments, given the 
strong tie between evaporation and sensible heat flux. Whereas the SWM experiments imposed 
the same diabatic heating anomaly in the different regions, Figure 5b shows that the effective 


strengths of these anomalies in the AGCM experiments varied across the country, with low 
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values in the western US and northeastern US (e.g., in the neighborhood of the Great Lakes) and 


with higher values toward the southeast. 


(iii) Basic Results. Figure 6 shows the responses of the 250-mb streamfunction, precipitation, 
and 2-m air temperature fields to the imposed dryness in Region S in Figure 2. (That is, the 
figure shows the difference between the mean fields for Experiment AGCM-S and the 
corresponding mean fields for the control.) In the figure, the solid lines, dashed lines, and heavy 
white lines correspond respectively to differences that are significantly different from zero at the 
95%, 99%, and 99.9% confidence levels, as determined by a t-test. Clearly seen in Figure 6a is 
an induced wave pattern in the streamfunction difference field, with a high centered just to the 
north of the dry area, a low further to the northeast, and another high even further to the 
northeast, covering the south of Greenland. This wave pattern is clearly evident in the 2-m 
temperature difference field (Figure 6c) and is also suggested in the precipitation difference field 
(Figure 6b). In the particularly affected areas to the north and east of Region S, precipitation is 
reduced by as much as 0.4 mm/day, and temperature is increased by as much as 1K — large 
differences given the 2-month averaging period. It is worth emphasizing here that, given the 
design of the experiment, the responses seen in Figure 6 stem solely from localized soil moisture 
drying in the AGCM - the imposed soil moisture anomalies are indeed having an impact on the 


atmospheric circulation and on remote near-surface meteorology. 


Results for all experiments are summarized in Figures 7, 8, and 9. Figure 7 shows the 
full streamfunction anomalies induced by imposed dryness in the different areas. The anomaly 
patterns naturally vary among the experiments, far more than they do among the corresponding 


SWM experiments. This is fully expected given the relative complexity of the AGCM. The 


17 


384 


385 


386 


388 


389 


392 


393 


395 


396 


397 


398 


400 


401 


405 


406 


differences between the SWM and AGCM experiments in their atmospheric circulation 
anomalies can be traced to their differences in stationary wave forcing anomalies. Relative to our 
SWM experiments, which utilize simple, localized diabatic heating anomalies (Figure 1), the 
AGCM experiments produce more realistic diabatic heating anomalies that vary regionally, and 
they also (unlike the particular SWM experiments we performed) account for anomalies in 
transient flux convergences. In addition, unlike the deterministic SWM results, the AGCM 
results are affected by chaotic atmospheric dynamics, necessitating an averaging over a large 
number of ensemble members, and even with this averaging, sampling error still has some 
impact on the AGCM results. Even so, the streamfunction anomalies produced in corresponding 
AGCM and SWM experiments do agree, at least to first order. Notice, for example, that 
experiments SWM-F and AGCM-F both produce a similar swath of high values across southern 
Canada, and experiments SWM-O and AGCM-O both show that dryness on the far east coast 
locates the high positive streamfunction values relatively far to the east. While the patterns 
produced in experiments SWM-P and AGCM-P and in SWM-Q and AGCM-Q differ 
significantly, this difference can be traced to the fact that in the AGCM, the zeroing of the 
rainfall in these regions did not lead to a large reduction in surface latent heat flux (and thus to a 
large increase in the heating of the overlying air through enhanced surface sensible heat flux), 
since soil moistures in the control simulations were already dry in these regions during the 
months considered (see Figure 5b). The SWM and AGCM experiments were thus fundamentally 


different for both of these regions. 


Consider now the experiments for which dryness is imposed within the area demarcated 
in Figure 4a (AGCM-B,C,D,K,L,M,R,S,T; the associated experimental results are enclosed 


within the red solid line in Figure 7). An interesting feature seen for all of these experiments is 
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the production of a high (with experiment-dependent magnitude) over the western-central US 
and the corresponding production of a low that straddles the northeastern coast. Though the 
details differ (again, the AGCM is much more complex than the SWM, and the AGCM results 
shown here are affected somewhat by sampling), the basic first-order pattern produced by an 
imposed land surface dryness in the interior region is quite similar to that found for the SWM. 
(This will be demonstrated further in section 3c below; averaging the precipitation changes 
produced by the AGCM experiments that focus on the demarcated interior region produces a 
field very similar to that shown in Figure 4b for the SWM experiments.) In other words, the 
AGCM results support, certainly to first order, the idea that dry conditions anywhere within the 


demarcated interior region produce a common atmospheric circulation response. 


Beyond providing information on atmospheric circulation, the AGCM experiments are 
particularly interesting because they provide information on precipitation and near-surface air 
temperature response to the imposed dryness anomalies — they provide, in particular, some 
indication of how surface dryness can affect remote near-surface meteorology. Such remote 
impacts are examined further in section 3d. For now, consider Figure 8, which shows the 
precipitation anomalies produced in the different AGCM experiments. To some degree, several 
of the experiments in the demarcated interior continental region produce similar precipitation 
reductions in southern Canada and in the eastern half of the US, down to the Gulf of Mexico 
(with details, of course, differing with experiment). Some similarities in precipitation response 
would indeed be expected from the aforementioned similar response in the streamfunction field, 
given that specific changes in large-scale atmospheric circulation can induce specific changes in 


precipitation. 
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Figure 9 shows the 2-meter air temperature anomalies produced in the AGCM 
experiments. With a few exceptions (namely, I, J, P, and Q — desert regions, for which the 
zeroing of precipitation did not strongly affect the soil moisture), the impact of the imposed 
drying is locally very strong, as might be expected given the associated reduced level of local 
evaporative cooling (Figure 5b). Many experiments, however, show in addition a remote air 
temperature impact. Several of the experiments focusing on dryness in the demarcated region in 
Figure 4a show warming along a swath extending from north-central North America to the Gulf 
of Mexico, a similarity that may stem from the aforementioned similarity in the streamfunction 


response. 


c. Signals Inherent in Observations-Based Data 


The results in section 3b relate to the AGCM’s inherent climate. While the hope is that 
they also represent nature, this is, of course, far from guaranteed. Unfortunately, verifying 
definitively our results with observations is made impossible by the limited extent of available 
observations-based data. To a large extent the shifts examined in the experiments above are 
subtle; given the large number of simulations constituting our ensembles, we are effectively 
examining here the net impact of shifts in probability distribution functions (PDFs). Reanalyses, 
which typically cover less than 40 years, are inadequate for a proper analysis of PDF shifts. 
Another key difficulty is the fact that when soil moisture in nature is dry in, say, region K of 
Figure 2, it also tends to be anomalously dry elsewhere (say, in region L and M), complicating 


tremendously the isolation of specific impacts associated with region K’s dryness. 
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The analyses above, however, suggest a way of addressing at least the second difficulty. 
Despite the spatial shifts in the local diabatic heating induced by the different dryness regions in 
Figure 2, the atmosphere in the different SWM and AGCM experiments responds in a similar 
way — dry conditions (or diabatic atmospheric heating associated with dry conditions) in the 
interior of the US tend to produce a similar streamfunction anomaly pattern. A possible 
inference from this result is that in the complex set of available observations, the key feature to 
look for is not a soil moisture deficit in a specific region but an average dryness in a large 
continental region (namely, that delimited in Figure 4a), without regard to the spatial distribution 
of this dryness — where it is maximized, where it is small, etc. If the SWM and AGCM results 
are valid, then the atmosphere should tend to show the same response to overall continental 


dryness regardless of the specific spatial pattern. 


With this in mind, we now analyze observations-based 250-mb stream function, 
precipitation, and air temperature fields associated with continental-scale antecedent land surface 


dryness. 


(i) Compositing strategy. We use a compositing approach to draw out signals, if they exist, in 
the observations-based data. Within each year of 1980-2014, we consider six 10-day periods: 
June 1-10, June 11-20, June 21-30, July 1-10, July 11-20, and July 21-30, amounting to a total of 
210 10-day periods. We tie to each period an antecedent soil moisture index in the following 
way. First, we standardize, on a daily basis, the MERRA-2 root zone soil moisture within each 
grid cell of the interior region demarcated in Figure 4a. Then, for the 10-day period in question, 
we identify the date falling fifteen days prior to the start of the period and compute, for that day, 


the average of the standardized soil moisture values across the region. A large negative value for 
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this average implies that, on average, the large-scale continental region was generally very dry 
15 days prior to the period in question. (See section 2c for a discussion of the realism of 
MERRA-? soil moisture values.) We use soil moistures at a 15-day lead rather than concurrent 
soil moistures to help ensure that the soil moisture signals on which we base our composites are 
not simply a passive reflection of the meteorology we are examining. Fifteen days lies above the 
typical time scales of planetary wave development and maintenance in response to local heating 


anomalies. 


Our composites are based on the driest 10% of antecedent average soil moisture values. 
From the 210 June-July decads examined, we thus compute the average 250-mb stream function, 


precipitation, and air temperature anomaly fields over the 21 driest cases. 


(ii) Results. Figure 10a, 10b, and 10c show respectively the resulting composites of 250-mb 
streamfunction (from ERA-Interim), precipitation (from CPCU, via MERRA-2), and 2-m air 
temperature (from ERA-Interim). The composited full streamfunction field shows a high along 
the U.S.-Canada border in the central part of the continent. The composited precipitation field 
shows a reduction of precipitation along the eastern border of the demarcated interior region 
along with some increases in precipitation in various surrounding areas, and the composited 


temperature field shows a positive temperature anomaly that spans much of the interior region. 


For comparison, Figure 10d shows the corresponding streamfunction anomaly field from 
the SWM experiments — it shows the sum of the fields obtained in the experiments imposing 
diabatic heating anomalies within the indicated region. The positioning of the positive 


streamfunction lobe from the observational composite agrees to first order with that in the SWM 
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results. While the observational composite is largely missing the negative lobe in the east, the 
square of the spatial correlation coefficient between the patterns in Figures 10a and 10d over the 


full area of each panel is reasonably high (1’=0.36). 


Finally, Figures 10e, 10f, and 10g show the corresponding fields derived from the AGCM 
experiments. To construct the figures, we simply add together the June-July anomalies of the 
experiments corresponding to the areas delimited by the heavy black lines, making the 
assumption (as we did above for the SWM) that, in the AGCM, the effects of the different 
dryness areas on continental meteorology combine linearly. (This assumption is strongly 
supported by the results of Koster et al. [2014]; see their Figure 5.) In comparing the 
observational composites and the AGCM results, the first thing to note is the difference in scale. 
This said, notice that the patterns inherent in the observational composites and the AGCM 
results are again similar. Only slight spatial shifts, for example, are seen in the positions of the 


large precipitation and air temperature anomalies. 


We emphasize here that these results do not prove that the SWM and AGCM experiments 
properly represent nature. Again, the experiments address the net impact of shifts in the 
underlying PDFs of the meteorological fields, shifts that cannot be captured by the limited 
observational record. That is, even if the model results represented nature perfectly (which, 
given their inherent biases, they presumably do not), we would still expect to see differences 
between the observational and model results. Nevertheless, the comparisons in Figure 10 are 
promising; they show that the experimental results are at least consistent with available 


observations. 
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d. Additional Findings Regarding Remote Land Surface Impacts 


The similarity in the pattern of the atmospheric response to dryness in different parts of 
the North American continent is the main result of this study. Nevertheless, the AGCM 
experiments hold a wealth of information that can be tapped for further understanding of land- 
atmosphere interactions. Here, we distill the extensive amount of information in Figures 7-9 into 
two simple aggregate quantities, quantities that together characterize the first-order impact of soil 
moisture anomalies on meteorological variables in remote locations and how this impact varies 


with location. 


The first aggregate quantity, the “impact extent”, measures the spatial extent of the 
impact of each imposed dryness region in Figure 2. While it is related in concept to more 
sophisticated diagnostics described in the literature for characterizing the remote reach of 
localized anomalies (e.g., the climate sensitivity maps constructed by Barsugli et al. (2006) from 
the results of SST warm patch experiments), the diagnostic we use here is more limited, focusing 
only on the spatial extent of this reach rather than on the amplitude of the response. Our quantity 
is computed as follows. For a given experiment, and for a given meteorological variable 
(precipitation or air temperature), we count the number of 1°x1° grid cells in North America (at 
least the part of North America shown in Figure 11) for which the change in the meteorological 
variable determined for the experiment was significantly different from zero at the 99% 
confidence level or higher, according to a t-test. This number, which indeed reflects the area 
within the dashed contours in Figures 8 and 9, is then plotted in Figure 11 — the number 
computed for a given experiment is plotted within the corresponding imposed dryness area for 


that experiment. 
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Figure 11a shows that for precipitation, the impact extent is highest for the dryness 
regions encompassing the Gulf Coast and parts of the Mississippi Valley (regions M,N, S, T). 
When these areas are forced to be dry, the remote (plus local) impacts of the dryness are felt 
across a relatively large area. High values also appear in an area around Montana (regions B and 
C). Values along either coast are relatively small, meaning that when these areas are dry, the 
spatial extent of the resulting impact on precipitation is slight. In fact, as might be expected, the 
impact extent varies to some degree with the impact of the drying on the June-July evaporation 
rates (Figure 5b). (Note that a 99% confidence level suggests that roughly 1%, or 36, of the 
3621 grid cells considered here will be falsely labeled significant. Values below 48 in the plots 


are shaded gray, suggesting little or no impact.) 


For air temperature (Figure 11b), region S has the greatest impact extent, followed by 
region B. Aside from region O, dryness regions along either coast again have minimal or no 
impact. Curiously, for some experiments, spatial impacts on air temperature appear smaller than 
those on precipitation, despite the fact that for air temperature, a local impact at the very least is 


almost guaranteed due to a reduction in local evaporative cooling. 


The second aggregate quantity, the “sensitivity to remote dryness”, measures the degree 
to which a given grid cell is affected by soil moisture anomalies in different regions. This 
sensitivity is computed at each grid cell as the number of different experiments for which a 
meteorological quantity (precipitation or air temperature) at the grid cell was modified at a 
significance level of 99% or higher, according to a t-test. Figure 12a shows the spatial 
distribution of the sensitivity for precipitation. Notice that 4 to 6 distinct dryness areas are able 
to affect the June-July precipitation in various grid cells along the eastern half of the US-Canada 


border; by this measure, these grid cells are particularly susceptible to remote effects. In general, 
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the eastern half of the US is more sensitive to US soil moisture dryness than the western half. 
The corresponding plot for air temperature (Figure 12b) also shows an east-west contrast but 
locates the region of highest “sensitivity to remote dryness” along the southeastern coast. Of 
course, it is important to keep in mind that the patterns in Figure 12 are determined, to some 
degree, by the distribution of dryness regions in Figure 2. If we had examined additional 
experiments with imposed dryness in Canada and Mexico, the patterns in Figure 12 would 


presumably be modified, particularly in or near these outside regions. 


The quantities plotted in Figures 11 and 12 are, in some ways, converses of each other; 
Figure 11 in essence shows the strength (local plus remote) of a “source” of meteorological 
anomalies, and Figure 12 shows, in a way, the strength of their “sink”. The locations of the 
sources and sinks are different, underlining the geographical connections that exist between land 


areas through the overlying circulation and its sensitivity to surface processes. 


4. Summary and Discussion 


The simulations with the stationary wave model (SWM) suggest that when a diabatic 
heating anomaly is imposed in the interior of the continental US (anywhere within the delimited 
area in Figure 4a) during boreal summer, the atmospheric circulation responds in a very similar 
way: it produces a positive eddy streamfunction anomaly in the west-central part of the 
continent and a negative anomaly to the east. This basic response is also seen, at least to first 
order, in AGCM experiments when soil moisture contents in different parts of the continental 
interior are forced to be dry. The AGCM experiments further allow us to identify impacts of the 


imposed soil moisture dryness on remote near-surface meteorological fields. As indicated by the 
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averages in the right-hand column of Figure 10, dry soil moistures in the interior continental 
region tend to produce warm conditions within the region (Figure 10g) and reduced rainfall a 
little to the east (Figure 10f). Observational composites (Figures 10a-c) are consistent, again to 
first order, with the modeled impacts. (It is worth emphasizing again that all of these results are 
for June-July only; land dryness impacts may be reduced during non-summer months due to 


reductions in surface turbulent fluxes associated with reduced solar forcing.) 


The overall behavior established in our experiments is strongly suggestive of a positive 
feedback loop, as captured in the schematic in Figure 13: dry conditions within the continental 
interior lead to changes in the atmospheric circulation that in turn lead to further warming and 
drying there. Consider, for example, the results of imposing soil moisture dryness in region L of 
Figure 1. According to Figure 71, this leads to a 250-mb streamfunction high over region L that 
also extends further to the north. Figures 81 and 91 show that this in turn leads to a warm and dry 
anomaly just north of the region — conditions that, according to experiment AGCM-D, should 
lead to additional warming and drying in the area (Figures 8d and 9d). The positive feedback 
loop outlined in Figure 13 may have impacts on the spatial extent of drought, helping a local dry 
anomaly in, say, Kansas and Nebraska spread spatially across the continental interior. The 


mechanics of this feedback loop are worthy of further study. 


What are the mechanisms by which an atmospheric circulation anomaly translates to 
anomalies in surface air temperature and precipitation? At these latitudes, streamfunction 
anomalies and geopotential height anomalies are largely coincident, and as discussed in Koster et 
al. (2014), an upper-level positive height anomaly can affect surface temperature through an 
associated increase in subsidence, which promotes cloudless skies and thereby an increased 


amount of incoming solar radiation reaching the surface. Lower level circulation changes, 
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however, are undoubtedly also important, as they can lead to precipitation anomalies through 
changing atmospheric moisture transport and convergence in the lower troposphere. We should 
note that the lower level impacts amongst the experiments are less systematic. Anomalies in 
850mb wind circulation patterns (not shown), as generated by the different AGCM experiments, 
do not show the level of agreement seen in Figure 7 largely because the anomalies are of smaller 
spatial extent and tend to be more centered over the individual imposed dryness regions. This 
said, however, some agreement in the 850mb wind response is seen, and supplemental analyses 
of the AGCM results indicate that the west-central North America streamfunction anomaly 
produced in most of our experiments does tend to tilt westward with height. This westward tilt 
may explain why the precipitation anomalies produced in our AGCM experiments tend to lie 
slightly eastward of the air temperature anomalies (compare Figures 10f,g), given that 
precipitation deficits are largely maintained by reduced atmospheric moisture flux convergences 
due to changes in low-level atmospheric circulation (e.g., Wang et al 2010; Wang and Schubert 
2014). Meanwhile, changes in low-level winds bring, in some cases, warm air in from the south, 
further increasing temperature (not shown). A full understanding of the mechanisms linking the 
3-D circulation patterns to surface anomalies needs additional study through focused 
experimentation and associated budget analysis for surface energy, air temperature and 


atmospheric moisture. 


While the feedback loop outlined in Figure 13 is intriguing, we must emphasize again a 
critical aspect of our results, namely, the need for large ensembles (96-192 members) to draw out 
the indicated patterns in Figures 7-9. The size of these ensembles speaks to the weakness of the 
model atmosphere’s response to the imposed surface anomalies; as noted above, our experiments 


effectively address shifts in the PDFs of atmospheric response. Dry conditions for a given year 
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in the central US will not necessarily lead to the indicated streamfunction anomalies during that 
year and an associated strengthening of the dry anomaly; all we can claim is an increased 


likelihood for this to happen. 


In other words, the feedback’s operation in the AGCM is probabilistic. It is presumably 
probabilistic in nature as well. The fact that the magnitudes of the responses in Figure 10 are 
larger for the observational composites than for the model results (while the patterns are roughly 
the same) in fact may suggest that the feedback in the model is weaker than it is in nature. 
Additional work is needed to learn if this is true and whether the feedback is indeed large enough 


in nature to be useful for, say, forecast applications. 


Perhaps most intriguing of all is our finding, through both the SWM and AGCM 
analyses, that the atmospheric circulation responds in such a consistent way to surface heating 
anomalies regardless of the specific geographical positioning of these anomalies within the 
continental interior (Figures 3 and 7). We are currently performing analyses aimed at explaining 
this marked similarity in response and the relatively reduced response produced when the dry 


anomalies are near the western or eastern coast. 
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Figure 1. a. Vertical profile of idealized diabatic heating anomaly imposed in the SWM 
atmosphere over a selected geographical area. b. Spatial distribution of imposed 
idealized heating anomaly near the ground surface (o = 0.9966), for experiment SWM-L. 
c. Vertical and zonal structure of the ensemble mean June-July diabatic heating 
anomalies produced in Experiment AGCM-L, averaged between 37°N and 44°N. Units: 


K/day. 


Figure 2. Areas of imposed dryness considered in the 21 SWM and AGCM experiments (labeled 
A through U). In a given AGCM experiment, the indicated region was forced to be dry 
through the artificial zeroing of the incident precipitation in April through June. The 
number of AGCM ensemble members associated with each experiment is indicated in 


parentheses. 


Figure 3. Eddy streamfunction anomaly fields at 250 mb produced in the SWM experiments. 


The boxes indicate the locations for which the idealized diabatic heating source is placed. 


Figure 4. Indication of the similarity between the 250-mb eddy streamfunction responses in the 
different SWM experiments. a. Location of the 9 experimental regions for which the 
square of the spatial correlation between the experiment’s 250-mb eddy streamfunction 
field and the average field over all 21 experiments over the North American area shown 
is 0.65 or higher (roughly 2/3 of variance explained by the averaged field). b. Average of 
the computed 250-mb eddy streamfunction fields over the nine experiments that examine 
heating within the region enclosed by heavy dotted lines in (a). c. Square of the spatial 


correlation (1?) between each experiment’s 250-mb eddy streamfunction field (from the 


37 


808 


809 


810 


811 


812 


813 


814 


815 


816 


817 


818 


820 


821 


823 


825 


826 


827 


829 


SWM analysis) and the 9-experiment average in (b). The value for a given experiment is 


plotted within the corresponding box in the figure. 


Figure 5. a. Percentile of average July soil moisture content in the different AGCM experiments, 


with rankings based on the full ensemble of control simulations. The percentile shown 
for a given grid cell is for the experiment in which that grid cell (along with its neighbors 
within the marked box) was forced to be dry during April through June; the map is thus a 
composite of results from different AGCM experiments. b. Difference (experiment 
minus control) in the June through July evaporation computed in the AGCM 
experiments. As in (a), the map shown is a composite; the difference shown for a given 
grid cell is for the experiment in which that grid cell (along with its neighbors within the 


marked box) was forced to be dry during April through June. 


Figure 6. a. Anomalies (10° m?/s) in June-July 250-mb streamfunction (experiment minus 


control) for Experiment AGCM-S. The imposed dryness in Region S has induced a 
wavelike pattern in the streamfunction field. b. Same, but for June-July precipitation 


anomalies (mm/day). c. Same, but for June-July 2-m air temperature anomalies (°K). 


Figure 7. a. Anomalies in June-July 250-mb streamfunction (experiment minus control) for 


Experiment A. Significance levels (according to a t-test) are shown as contours, with 
solid black for the 95% confidence level, dotted for the 99% confidence level, and solid 
white for the 99.9% confidence level. b-u: Same, but for Experiments B-U. The heavy 
red line encloses the experiments corresponding to the continental interior region 


demarcated in Figure 4a. 


Figure 8. Same as Figure 7, but for precipitation anomalies (mm/day). 
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Figure 9. Same as Figure 7, but for 2-m air temperature anomalies (°K). 


Figure 10. a. Composited field of 250-mb streamfunction anomalies from the ECMWF ERA- 


Interim reanalysis, built from the 10% of the 210 June and July decads during 1980-2014 
with the driest antecedent (15 day lead) continental soil moistures in the region indicated. 
b. Same, but for MERRA-2 rain-gauge corrected precipitation. c. Same, but for ERA- 
Interim 2-m air temperatures. d. Average of the June-July 250-mb eddy streamfunction 
anomalies produced in the 9 SWM experiments indicated. e. Sum of the June-July 250- 
mb streamfunction anomalies produced in the 9 AGCM experiments indicated. f. Same, 


but for precipitation. g. Same, but for 2-m air temperature. 


Figure 11. a. Number of 1°x1° grid cells across the indicated region (130W-60W, 20N-70N) for 


which precipitation is modified at the 99% confidence interval by the dryness imposed in 
each experiment. (The number of grid cells affected in a given experiment is plotted 
within the dryness region associated with that experiment.) Gray shading indicates very 


roughly the number expected by chance. b. Same, but for 2-m air temperature. 


Figure 12. a. For each 1°x1° grid cell, the number of dryness regions (as outlined in the heavy 


black lines) that, in our experiments, induce precipitation anomalies significant at the 


99% confidence level. b. Same, but for 2-m air temperature. 


Figure 13. Schematic of a potential feedback loop suggested by the AGCM results. The dry soil 


moistures induce (1) a specific circulation pattern that in turn induces (2) additional 
warming and drying over the central US. This in turn leads (3) to higher sensible heat 
fluxes from the land surface in the central US, which can enhance (4) the atmospheric 
circulation anomaly that originally produced the warming and drying. 
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Figure 1. a. Vertical profile of idealized diabatic heating anomaly imposed in the SWM 
atmosphere over a selected geographical area. b. Spatial distribution of imposed idealized 
heating anomaly near the ground surface (o = 0.9966), for experiment SWM-L. c. Vertical and 
zonal structure of the ensemble mean June-July diabatic heating anomalies produced in 
Experiment AGCM-L, averaged between 37°N and 44°N. Units: K/day. 
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Figure 2. Areas of imposed dryness considered in the 21 SWM and AGCM experiments (labeled 
A through U). In a given AGCM experiment, the indicated region was forced to be dry through 
the artificial zeroing of the incident precipitation in April through June. The number of AGCM 
ensemble members associated with each experiment is indicated in parentheses. 
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SWM Results: Eddy Streamfunction (10° m’/s) 


Figure 3. Eddy streamfunction anomaly fields at 250 mb produced in the SWM experiments. 
The boxes indicate the locations for which the idealized diabatic heating source is placed. 
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Figure 4. Indication of the similarity between the 250-mb eddy streamfunction responses in the 
different SWM experiments. a. Location of the 9 experimental regions for which the square of 
the spatial correlation between the experiment’s 250-mb eddy streamfunction field and the 
average field over all 21 experiments over the North American area shown is 0.65 or higher 
(roughly 2/3 of variance explained by the averaged field). b. Average of the computed 250-mb 
eddy streamfunction fields (10° m7/s) over the nine experiments that examine heating within the 
region enclosed by heavy dotted lines in (a). c. Square of the spatial correlation (1?) between 
each experiment’s 250-mb eddy streamfunction field (from the SWM analysis) and the 9- 
experiment average in (b). The value for a given experiment is plotted within the corresponding 
box in the figure. 
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Figure 5. a. Percentile of average July soil moisture content in the different AGCM experiments, 
with rankings based on the full ensemble of control simulations. The percentile shown for a 
given grid cell is for the experiment in which that grid cell (along with its neighbors within the 
marked box) was forced to be dry during April through June; the map is thus a composite of 
results from different AGCM experiments. b. Difference (experiment minus control) in the June 
through July evaporation computed in the AGCM experiments. As in (a), the map shown is a 
composite; the difference shown for a given grid cell is for the experiment in which that grid cell 
(along with its neighbors within the marked box) was forced to be dry during April through June. 
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Figure 6. a. Anomalies (10° m?/s) in June-July 250-mb streamfunction (experiment minus 


control) for Experiment AGCM-S. Significance levels (according to a t-test) are shown as 
contours, with solid black for the 95% confidence level, dotted for the 99% confidence level, and 
solid white for the 99.9% confidence level. The imposed dryness in Region S has induced a 
wavelike pattern in the streamfunction field. b. Same, but for June-July precipitation anomalies 
(mm/day). c. Same, but for June-July 2-m air temperature anomalies (°K). 
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Figure 7. a. Anomalies in June-July 250-mb streamfunction (experiment minus control) for 
Experiment A. Significance levels (according to a t-test) are shown as contours, with solid black 
for the 95% confidence level, dotted for the 99% confidence level, and solid white for the 99.9% 
confidence level. b-u: Same, but for Experiments B-U. The heavy red line encloses the 
experiments corresponding to the continental interior region demarcated in Figure 4a. 
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Figure 8. Same as Figure 7, but for precipitation anomalies (mm/day). 
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Figure 9. Same as Figure 7, but for 2-m air temperature anomalies (°K). 
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Figure 10. a. Composited field of 250-mb streamfunction anomalies from the ECMWF ERA- 
Interim reanalysis, built from the 10% of the 210 June and July decads during 1980-2014 with 
the driest antecedent (15 day lead) continental soil moistures in the region indicated. b. Same, 
but for MERRA-2 rain-gauge corrected precipitation. c. Same, but for ERA-Interim 2-m air 
temperatures. d. Average of the June-July 250-mb eddy streamfunction anomalies produced in 
the 9 SWM experiments indicated. e. Sum of the June-July 250-mb streamfunction anomalies 
produced in the 9 AGCM experiments indicated. f. Same, but for precipitation. g. Same, but for 
2-m air temperature. 
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Figure 11. a. Number of 1°x1° grid cells across the indicated region (130W-60W, 20N-70N, or 
3621 cells) for which precipitation is modified at the 99% confidence interval by the dryness 
imposed in each experiment. (The number of grid cells affected in a given experiment is plotted 
within the dryness region associated with that experiment.) Gray shading indicates very roughly 
the number expected by chance. b. Same, but for 2-m air temperature. 
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Figure 12. a. For each 1°x1° grid cell, the number of dryness regions (as outlined in the heavy 
black lines) that, in our experiments, induce precipitation anomalies significant at the 99% 
confidence level. b. Same, but for 2-m air temperature. 


120W 


60W 


120W 


51 


60W 


Regional soil moisture anomalies 


imposed in US interior 


1 


Atmospheric circulation anomalies 
A high over west-central US and a low to the east 


i) 


4 


Enhance/adjust 


Surface meteorological (T2m and P) anomalies: 
warming and dryness over central US 


2 


Diabatic heating anomalies: 


Additional surface sensible heating over central US 


Figure 13. Schematic of a potential feedback loop suggested by the AGCM results. The dry soil 
moistures induce (1) a specific circulation pattern that in turn induces (2) additional warming and 
drying over the central US. This in turn leads (3) to higher sensible heat fluxes from the land 
surface in the central US, which can enhance (4) the atmospheric circulation anomaly that 
originally produced the warming and drying. 
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